Relationship between depth-wise refractive index and biomechanical properties of human articular cartilage

Abstract. Significance Optical properties of biological tissues, such as refractive index (RI), are fundamental properties, intrinsically linked to the tissue's composition and structure. We hypothesize that, as the RI and the functional properties of articular cartilage (AC) are dependent on the tissue’s structure and composition, the RI of AC is related to its biomechanical properties. Aim This study aims to investigate the relationship between RI of human AC and its biomechanical properties. Approach Human cartilage samples (n=22) were extracted from the right knee joint of three cadaver donors (one female, aged 47 years, and two males, aged 64 and 68 years) obtained from a commercial biobank (Science Care, Phoenix, Arizona, United States). The samples were initially subjected to mechanical indentation testing to determine elastic [equilibrium modulus (EM) and instantaneous modulus (IM)] and dynamic [dynamic modulus (DM)] viscoelastic properties. An Abbemat 3200 automatic one-wavelength refractometer operating at 600 nm was used to measure the RI of the extracted sections. Similarly, Spearman’s and Pearson’s correlation coefficients were employed for non-normal and normal datasets, respectively, to determine the correlation between the depth-wise RI and biomechanical properties of the cartilage samples as a function of the collagen fibril orientation. Results A positive correlation with statistically significant relations (p−values<0.05) was observed between the RI and the biomechanical properties (EM, IM, and DM) along the tissue depth for each zone, e.g., superficial, middle, and deep zones. Likewise, a lower positive correlation with statistically significant relations (p−values<0.05) was also observed for collagen fibril orientation of all zones with the biomechanical properties. Conclusions The results indicate that, although the RI exhibits different levels of correlation with different biomechanical properties, the relationship varies as a function of the tissue depth. This knowledge paves the way for optically monitoring changes in AC biomechanical properties nondestructively via changes in the RI. Thus, the RI could be a potential biomarker for assessing the mechanical competency of AC, particularly in degenerative diseases, such as osteoarthritis.

Aim: This study aims to investigate the relationship between RI of human AC and its biomechanical properties.
Approach: Human cartilage samples (n ¼ 22) were extracted from the right knee joint of three cadaver donors (one female, aged 47 years, and two males, aged 64 and 68 years) obtained from a commercial biobank (Science Care, Phoenix, Arizona, United States).The samples were initially subjected to mechanical indentation testing to determine elastic [equilibrium modulus (EM) and instantaneous modulus (IM)] and dynamic [dynamic modulus (DM)] viscoelastic properties.An Abbemat 3200 automatic one-wavelength refractometer operating at 600 nm was used to measure the RI of the extracted sections.Similarly, Spearman's and Pearson's correlation coefficients were employed for non-normal and normal datasets, respectively, to determine the correlation between the depth-wise RI and biomechanical properties of the cartilage samples as a function of the collagen fibril orientation.
Results: A positive correlation with statistically significant relations (p − values < 0.05) was observed between the RI and the biomechanical properties (EM, IM, and DM) along the tissue depth for each zone, e.g., superficial, middle, and deep zones.Likewise, a lower positive correlation with statistically significant relations (p − values < 0.05) was also observed for collagen fibril orientation of all zones with the biomechanical properties.

Conclusions:
The results indicate that, although the RI exhibits different levels of correlation with different biomechanical properties, the relationship varies as a function of the tissue depth.This knowledge paves the way for optically monitoring changes in AC biomechanical properties nondestructively via changes in the RI.Thus, the RI could be a potential biomarker for assessing the mechanical competency of AC, particularly in degenerative diseases, such as osteoarthritis.

Introduction
Articular cartilage (AC) is a specialized connective tissue covering the ends of bones in diarthrodial joints, such as the knee, hip, and shoulders. 1The primary function of AC is the provision of frictionless articulation and transmission of physical loading to the underlying subchondral bone without leading to high stresses. 2Mature AC has limited repair capability due to the lack of vasculature and nerves; thus, it is susceptible to degenerative changes driven by age or injury.The extracellular matrix (ECM) of AC mainly comprises ∼5% to 10% proteoglycan (PG) macromolecules, 10% to 20% collagen fibril meshwork, and around 60% to 80% water. 3he structure and composition of AC vary along the tissue depth. 4,5Structurally, AC can be characterized into four distinctive zones based on collagen fibril orientation.These zones are the superficial zone (SZ), middle zone (MZ), deep zone (DZ,) and calcified zone (CZ).The SZ comprises ∼0% to 10% of AC overall thickness, 6 with collagen fibers aligned parallel to the articular surface of the tissue.The SZ is in direct contact with synovial fluid present in the joint space and protects the ECM from shear stresses.The transitional zone or MZ covers ∼10% to 30% of AC thickness and plays a key role in load transmission.Collagen fibers in the MZ exhibit random orientation.The DZ covers ∼70% to 90% of the ECM 6 and interfaces with CZ.The DZ provides the highest resistance to compressive loading 2 and has collagen fibers oriented perpendicular to the cartilage-bone interface.The water content decreases, the PG content increases, and the density of the collagen fibril network increases from the SZ to the DZ. 5,7,8urthermore, in the SZ, chondrocytes are flattened, whereas they exhibit a spherical shape in the MZ and DZ.In the DZ, chondrocytes have a columnar cluster arrangement with vertical orientation, parallel to the collagen fibers. 5The SZ holds a relatively higher number of chondrocytes compared with the MZ and DZ. 5 Compositionally, PG macromolecules are responsible for cartilage compressive properties, whereas collagen fibrils contribute more to the dynamic and tensile properties. 9,10Alteration of the depth-wise profile of the AC matrix constituents often results in progressive degeneration, leading to osteoarthritis (OA), one of the most prevalent degenerative joint diseases, with a global prevalence estimated to be up to 28% in adults more than 60 years old. 11OA is a multifaceted joint disease primarily linked with disruption of collagen fibers network, 12 loss of PG content, increased matrix water content, and subsequent compromise of the mechanical function of the tissue. 2he refractive index (RI) is one of the primary optical properties that characterizes the manner in which light interacts with a material as it traverses through it. 13The RI of biological tissues is a key parameter that contributes to the understanding of light propagation through the tissue, potentially providing diagnostic information on tissue integrity 2 and providing insight into their molecular composition and structural integrity.Previous studies have highlighted the importance of polarized light interactions in characterizing these optical properties. 14ecent studies 13,15,16 have reported that both absorption and scattering of the incoming light, combined with the RI, play critical roles in understanding the interaction of light in biological tissues.For example, tumor development in biological tissues is often characterized by significant alteration in the morphology and density of cells in the pathological tissue; these structural changes impact the light scattering patterns, which can subsequently be used as a biomarker for tissue characterization (e.g., healthy versus cancerous). 15In addition, studies have presented key findings related to the variation of RI values across samples from human liver, kidney, myocardium, and skin, 16 suggesting the RI as a potential optical property for tissue characterization.The RI has shown significant potential for the characterization and diagnostic assessment of tissue pathologies. 13As AC is a highly anisotropic, optically turbid, and inhomogeneous tissue, 17 its scattering nature strongly influences the propagation of light, which is further determined by the tissue composition, integrity, and structure along the optical path. 18iomechanical properties of AC originate from the multi-phasic heterogeneous structure of the ECM.The constituents of the ECM and their interactions determine the non-linear mechanical response of AC. 19,20 PGs are the primary regulator of AC mechanical behavior under the scenario in which the interstitial fluid flow has seized. 21Furthermore, both the collagen fiber meshwork and the interstitial fluid of the ECM control the mechanical behavior of AC during cyclic and transient loads. 22OA is known to have significant impacts on the AC structure and function during the early disease stages.For example, PG loss is progressive and results in the reduction of cartilage load-bearing capacity in equilibrium. 23Subsequently, disease progression leads to surface fibrillation due to collagen network damage, and a portion of the tissue swellingrelated pretension of the fibers is compromised.5][26] However, scientific studies on the effect of alteration of human AC biomechanical properties on its optical properties are scarce.As far as we know, no study has investigated the relationship between the human AC RI and its biomechanical properties.
In this study, we hypothesize that there is a relationship between the depth-wise RI of human AC and its biomechanical properties.This hypothesis is based on our initial findings that the RI of human AC varies with its ECM structure and composition. 27Furthermore, the scientific basis for the hypothesis is supported by the depth-dependent variation of AC major constituents, e.g., collagen fibers, chondrocytes density and shape, water distribution, and PG content.Consequently, it is our position that the depth-dependent composition and structure of human AC will result in varying light-tissue interaction along the tissue depth, resulting in varying relationships between the biomechanical properties of the tissue and RI along the tissue depth.To test this hypothesis, the RIs of tissue sections extracted from different depths of human patellar articular cartilage samples were measured, and their relationship with the bulk tissue biomechanical properties was investigated.

Specimens Preparations
Human cartilage samples (n ¼ 22) were extracted from the right knee joint of three cadaver donors (one female, aged 47 years, and two males, aged 64 and 68 years) obtained from a commercial biobank (Science Care, Phoenix, Arizona, United States).Sample collection was approved by the Research Ethics Committee of the Northern Savo Hospital District (Kuopio University Hospital, Kuopio, Finland, #134//2015).The numbers of samples obtained from the three donors were eight, six, and eight.The samples were extracted from the central and anterior parts of the medial femur (n ¼ 5 and 3, respectively) and the central and anterior parts of the lateral tibia (n ¼ 6 and 4, respectively).The samples were extracted as cylindrical osteochondral plugs using a dental drill (diameter = 8 mm).The samples were initially subjected to mechanical indentation testing to determine the elastic [equilibrium modulus (EM) and instantaneous modulus (IM)] and dynamic [dynamic modulus (DM)] viscoelastic properties and then divided into two specimens: A and B as shown in Fig 1 .Specimen A was subjected to RI measurement, whereas histological assessment was performed on specimen B (already reported in our previous study 27 ).The thickness of the cartilage in specimen A was determined as the mean value of thickness measurement from the sides of the specimen, using an optical microscope (Zeiss, STEMI, SV8, Oberkochen, Germany). 28pecimen A was then immersed in phosphate-buffered saline (PBS) and stored at −20°C before undergoing the cryosectioning procedure.To measure the RI of different cartilage layers, tissue sections (200 μm) were extracted from depths corresponding to the different zones of the articular cartilage (SZ, MZ, and DZ) using a cryosection procedure.In summary, the bone end of each plug was flattened and adjusted to be parallel with the articular surface using a bandsaw.The flattened bone end was then attached to a disc using a fixative medium (optimal cutting temperature compound, Thermo Fisher Scientific Ltd., Runcorn, United Kingdom).The disc was then mounted in the cryostat chamber (Leica CM3050 S, Leica Biosystems, Wetzlar, Germany) for ∼5 min at −20°C to undergo rapid freezing, followed by placement on the microtome within the chamber.The microtome has a cutting range of 0.5 to 300 μm, with a precision orientation of tissue sections of ∼8 deg.Prior to sectioning, a disposable blade (MX35 Ultra, Thermo Fisher Scientific Ltd., Runcorn, United Kingdom) with a thickness of 0.24 mm and a cutting angle of 34 deg was set inside the knife holder of the cryostat chamber, and the first 50 μm of the articular surface was cut and discarded to create a uniformly flat surface for sectioning.The removal of this initial section of the articular surface is due to the naturally curvy surface of the tissue, thus ensuring the flatness of the subsequently extracted sections.After flattening the articular surface, 10 consecutive 200-μm-thick sections were extracted and mounted on glass slides (thickness = 1 mm, Menzel-Gläser frosted microscope slides, Thermo Fisher Scientific Oy, Vantaa, Finland) and then stored for RI measurement in a humid box at 4°C.In total, 220 sections were extracted for the analysis.It is important to mention that the biomechanical properties of articular cartilage are not significantly altered by the freeze-thaw cycle, as established in previous studies, including research conducted by our group. 29Therefore, performing mechanical measurements on fresh cartilage and subsequent optical measurements post-freezing is a valid and widely accepted approach.
Histological evaluation was performed on the B specimens.The samples were fixed and decalcified in a solution containing formaldehyde and ethylene diamine tetra acetic acid for 21 days, followed by processing and embedding in paraffin for sectioning. 30For the measurement of collagen fibril orientation, 5-μm histological sections were cut perpendicular to the articular surface.Likewise, the 3-μm histological sections were cut perpendicular to the articular surface and stained with Safranin O for digital densitometry and Mankin, and Osteoarthritis Research Society International (OARSI) scoring, being used in our previous study. 27dentification of cartilage zones was based on the depth-wise orientation of the collagen fibers derived from polarized light microscopy (PLM).Thus, the SZ, MZ, and DZ were defined as collagen fibril orientation of approximately 20 deg or less, >20 and <70 deg, and >70 deg, respectively.Based on the literature, 28,30,31 the SZ, MZ, and DZ correspond to 0% to 10%, 10% to 30%, and 30% to 100% of the overall cartilage thickness, respectively.During the preparation and extraction of the samples, the tissue surface was kept hydrated at regular intervals with PBS containing protease inhibitors. 32

Indentation Testing
To measure the mechanical properties of the cartilage, an indentation test with a uniform planended indenter (1 mm diameter) was performed on the osteochondral samples (8 mm diameter), before splitting the sample into two.EM was determined using stepwise indentation stress relaxation tests and calculated from the slope of the linear least-squares fit to the equilibrium stressstrain points.The sample was mounted on a rigid sample holder followed by compression with the plane-ended indenter, with a pre-stress of 12.5 kPa followed by a four-step stress relation test. 33Each step consisted of 5% strain and 15 min of relaxation time.Based on a previous study, which reported an optically measured Poisson's ratio for bovine tibial cartilage, 34 the Poisson's ratio at equilibrium was set to 0.3. 19Similarly, the IM of cartilage was calculated from each peak stress point.Thus, we can obtain data points for the IM as a function of the applied strain.Due to the high loading rate, the samples were assumed to be incompressible, and the Poisson's ratio was set to 0.57. 35The obtained strain-dependent instantaneous moduli were corrected using Hayes equation, 36 which takes into account the measurement geometry.The DM was determined from a dynamic sinusoidal test, conducted using 2% strain amplitude with frequencies of 0.005, 0.05, 0.1, 0.25, 0.5, 0.625, 0.833, and 1 Hz.The DM was defined as the ratio of stress and strain amplitudes in the dynamic sinusoidal test. 19During dynamic testing, the samples were assumed incompressible, having Poisson's of 0.5. 35

Refractive Index Measurement
An Abbemat 3200 automatic one-wavelength refractometer operating at 600 nm was used to measure the RI of the extracted sections from specimen A. The device is equipped with an automatic temperature controller, extending precise and fast temperature control of the section with an uncertainty of AE0.01 K. [37][38][39] Before measurement, the calibration of the refractometer was performed using double distilled water. 38After placing the section on the sample holder, a magnetic cover was placed on the prism to cover it entirely to ensure better light-tissue contact without external stray light.To assess the device sensitivity and precision, the RI measurement of each tissue section was repeated three times, and the average of these three measurements was used to determine the RI of each section.All measurements were conducted at room temperature (around 22°C).

Collagen Orientation Measurement
The sections extracted from specimen B (Fig. 1) were subjected to additional analysis to determine the orientation of the collagen fibers using PLM.PLM imaging was performed using an Abrio PLM system (CRi, Inc., Woburn, Massachusetts, United States) staged on a conventional light microscope (Nikon Diaphot TMD, Nikon, Inc., Shinagawa, Tokyo, Japan). 40,41The Abrio PLM system consisted of a circular polarizer, a green bandpass filter, and a computer-controlled analyzer consisting of two liquid crystal polarizers and a charge-coupled device camera.All samples were imaged in uniform orientation with a 4.0× objective, which resulted in a pixel size of 2.53 × 2.53 μm. 2 In the orientation images, 0 and 90 deg correspond to the orientations parallel and perpendicular to the cartilage surface, respectively. 42

Statistical Analysis
To assess the relationship between the biomechanical properties of human articular cartilage and its depth-wise RI, the extracted sections from all samples were grouped based on the zonal segmentation across the tissue depth, e.g., the SZ, MZ, and DZ.The normality of the RI datasets was also checked prior to analysis, using the Anderson-Darling test.The normality test was performed on each group of the RI, PLM, and biomechanical properties datasets, and based on the outcome, the non-parametric Kruskal-Wallis test was used for the non-normal dataset.Similarly, Spearman's and Pearson's correlation coefficients were employed for nonnormal and normal datasets, respectively, to determine the correlation between the depth-wise RI and biomechanical properties of the cartilage samples as a function of collagen fibril orientation.
In the statistical tests, a p − value lower than 0.05 was considered to reflect statistical significance.All statistical analyses were performed using MATLAB (MathWorks, Natick, Massachusetts, United States, version R2019b).
A positive correlation with statistically significant relations (p − values < 0.05) was observed between the RI and the biomechanical properties (EM, IM, and DM) along the tissue depth for each zone (Table 1).Likewise, a lower positive correlation with statistically significant relations (p − values < 0.05) was also observed for collagen fibril orientation of all zones with the biomechanical properties (Table 2).The empty cells for Spearman's and Pearson's correlation reflect that the datasets of the variables concerned are non-normal and normal, respectively.In addition, the relationship between the zonal RI and biomechanical properties (IM, EM, and DM) shows an increasing trend for the SZ, MZ, and DZ.A similar trend can be observed with the collagen fibril orientation as illustrated in Fig. 3. Statistically significant differences were observed between the RI of the SZ and EM, IM, and DM, with Spearman's correlation of 0.4385, 0.2796, and 0.3003, respectively, and p < 0.001 in all cases, as illustrated in Table 1.Similarly, statistically significant differences were observed between the RI of the MZ and EM, IM, and DM, with Pearson's coefficient of 0.0615, 0.5107, and 0.4803, respectively, and p < 0.001 in all cases (Table 1).Likewise, statistically significant differences were noticed between the RI of the DZ and EM, IM, and DM, with Spearman's correlation of 0.3882 for EM and Pearson coefficient of 0.6152 and 0.5158 for IM and DM, respectively, and p < 0.001 in all cases as illustrated in Table 1.
Similar to RI, statistically significant correlations were observed between collagen fibril orientation in the SZ and EM, IM, and DM, with Spearman's correlation of 0.6738 for EM and Pearson coefficient of −0.3063 and 0.0938 for IM and DM, respectively, with p < 0.001 in all cases (Table 2).Similarly, statistically significant correlations were also noticed between collagen fibril orientation in the MZ and EM, IM, and DM, with Spearman's correlation of −0.0802, −0.2250, and 0.0506, respectively, and p < 0.001 in all cases (Table 2).Finally, statistically significant correlations were observed between collagen fibril orientation in the DZ with EM, IM, and DM, having Spearman's correlation = 0.4205 with EM and Pearson coefficients = 0.5220, and 0.2729, well within the range reported, 43 with IM and DM, respectively, and p < 0.001 in all cases as illustrated in Table 2.

Discussion
In this study, we investigated the relationship between the depth-wise RI and mechanical properties of human AC.Although the potential of RI for the characterization of articular cartilage integrity has been previously demonstrated, 1,44 no study, to the best of our knowledge, has investigated the relationship between the depth-wise RI and the biomechanical properties of human AC.As the propagation of light in biological tissues is impacted by the ECM constituents 45 and cartilage structure and composition are known to vary along its depth, we hypothesized that there would be a depth-dependent relationship between the cartilage RI and its mechanical characteristics. 46he observed variations in the relationship between the RI and biomechanical properties of human AC are probably related to depth-wise variation in the tissue structure, composition, and integrity.For example, collagen fibril orientation in the SZ becomes less parallel to the tissue surface in the early stages of degeneration in comparison to healthy tissue. 47This potentially leads to alteration of the RI as a function of tissue depth and degeneration. 27he knowledge of the depth-wise variation of RI and in relation to tissue properties and degenerative state is important for developing models for understanding the propagation of light in biological tissues, such as Monte Carlo simulation-a gold standard method for investigating the interaction of light in biological tissue. 48Thus, detailed knowledge of the variation RI with respect to the structure and composition of articular cartilage along its depth could provide critical information for accurately modeling the interaction of light in the tissue.
The relationship between the RI and the biomechanical properties of AC is likely due to an indirect relationship with its collagen fiber network. 31The collagen fibril network is known to have a major effect on cartilage mechanical properties, especially under dynamic and instantaneous loading. 19,49Studies have shown that the variation in the fibrillar architecture of different cartilage zones affects its mechanical properties. 49,50Furthermore, given its amount (22% of ECM's wet weight), varying depth-wise orientation, and fibril diameter, collagen is the major light-scattering component of the articular cartilage. 51,52Although chondrocytes may also contribute to the scattering of light in the cartilage, their effect is likely to be minimal due to their dispersed nature and relatively low concentration (2% of volume) in ECM.
The relationship between the RI and biomechanical properties of human AC provides an indication of the potential of the RI for diagnostic characterization of cartilage properties.This relationship is valuable as it can be further developed into an approach for non-destructively estimating the biomechanical properties of articular cartilage, which are important biomarkers of articular cartilage degeneration.For example, cartilage integrity based on the OARSI score has been shown to be related to the biomechanical properties of the tissue. 53Likewise, an increase in the OARSI score has been reported to be related to the decrease in the elastic (equilibrium/ dynamic/IM) and increase in the viscous (phase difference) properties of the cartilage. 21urthermore, it has been shown that healthy cartilage has significantly higher EM compared with that of early and advanced OA cartilages. 20,54Thus, the RI could provide an indicator for estimating articular cartilage integrity through the relation between the tissue's mechanical properties and its integrity score, evaluated via OARSI scoring.
Indentation testing has shown a good capability to reveal alterations in tissue properties in the SZ of cartilage, where the early signs of OA can be observed 47,55 or even initiated during traumatic joint injury (post-traumatic OA).This is often characterized by fibrillation of the superficial collagen network (one of the earliest signs of OA 56 ), further supporting the use of an indentation test in this study.Likewise, earlier studies have demonstrated the loss of PGs as a potential indicator for assessing the development of OA. [57][58][59] A decrease in cartilage PG content decreases its swelling pressure, resulting in increased permeability and alteration of the load-bearing ability of its matrix. 60Thus, upon PG content depletion, the permeability of the ECM would increase significantly. 61,62Consequently, a recent study 63 has reported that changes in permeability result in altered mechanical properties, leading to cartilage matrix damage and eventually OA and alteration of collagen integrity of the tissue. 19hus, the relationship between RI and biomechanical properties of articular cartilage could potentially enable the development of optical diagnostic methods that could enable the detection of early degenerative changes in the cartilage. 64,65ikewise, in accordance with the structure-function relation, the RI and the biomechanical response of biological tissues are intrinsically related to the tissue's structure and composition and thus could be a potential biomarker of tissue integrity.For example, the RI could serve as a useful indicator in cartilage tissue engineering for monitoring the deposition of cartilage ECM during tissue growth in a bioreactor.There are currently no studies in the literature, to the best of our knowledge, on the relationship between depth-wise RI and biomechanical properties of human articular cartilage with respect to potential in vivo application.
Taken together, the estimation of cartilage integrity via the prediction of the biomechanical properties from the RI (i.e., based on the relationship between OARSI scores and biomechanical properties of articular cartilage) can be potential applications of this study in the future.The main limitation of this study is the relatively low number of samples extracted from the knees of three cadavers.This limited our ability to draw a firm conclusion; thus, more detailed future investigations with a greater number of samples are needed.Samples extracted from cadavers of varying ages and from joints of varying degenerative states would allow for a better assessment of the relationship between the depth-wise RI and mechanical properties of human AC.The other key limitation of this study was the removal of the initial first few tens of micrometers from the articular surface.This was done to flatten the naturally curved articular surface to allow the preparation of uniformly flat sections.We believe that this did not significantly affect the estimated RI of the SZ in the case of healthy human articular cartilage samples; however, this may not be the case for pathological samples as it may result in the loss of some useful details in the SZ.

Conclusion
In conclusion, the present results indicate that there is a correlation between the RI and depthwise variation of biomechanical properties of human articular cartilage.This finding is consistent with the depth-wise variation in cartilage composition and structure and provides critical information that could be useful for accurate modeling of light-tissue interaction using methods such as Monte Carlo simulation.For example, accounting for the variation in the RI for the different zones of articular cartilage could enable the development of more accurate models of light propagation in cartilage as compared with using a single RI for the bulk tissue.Moreover, the outcome of this study could support the development of complex Monte Carlo simulation models that can take into account the depth-wise variation of the RI for forecasting cartilage degradation and progression of OA based on a prediction of the biomechanical properties.

Fig. 1
Fig.1Schematic illustration of the experimental setup for sample preparation.Horizontal sections were extracted from specimen A for RI measurement, and vertical sections were extracted from specimen B for histological assessment and microscopy (PLM).Mechanical indentation testing was performed to measure elastic (EM and IM) and dynamic viscoelastic (DM) properties of the samples.

4. 1
Mapping of Zonal Refractive Index, Biomechanical Properties (EM, IM, and DM) with Collagen Fibril Orientation of Human Articular Cartilage Depth-wise distribution of the RI, collagen fibril orientation, and biomechanical properties (EM, IM, and DM) is illustrated in the histograms in Fig. 2. The collagen fibril orientation shows a consistent depth-dependent variation across the tissue depth as shown in Fig. 3.The collagen fibril orientation for the SZ, MZ, and DZ shows variation between 11.7 and 59.6 deg, 21.5 and 85.7 deg, and 54.4 and 83.7 deg, respectively.

Fig. 2
Fig. 2 Zonal histograms for the RI, collagen fibril orientation (PLM), and histogramic illustration for distribution of biomechanical properties (EM, IM, and DM) of human articular cartilage.

Fig. 3
Fig. 3 Mapping of the zonal RI and biomechanical properties (EM, IM, and DM) with collagen fibril orientation [PLM angle (deg)] of human articular cartilage.

Table 1
Summary of correlations between the RI and cartilage biomechanical properties in the SZs, MZs, and DZs of the articular cartilage.

Table 2
Summary of correlation between the collagen fibril orientation angle and cartilage biomechanical properties in the SZs, MZs, and DZs of articular cartilage.